B(E1) Strengths from Coulomb Excitation of n Be 
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The B(El;l/2 + -> 1/2") strength for n Be has been extracted from intermediate energy Coulomb 
excitation measurements, over a range of beam energies using a new reaction model, the extended 
continuum discretized coupled channels (XCDCC) method. In addition, a measurement of the 
excitation cross section for 11 Be+ 208 Pb at 38.6 MeV/nucleon is reported. The _B(E1) strength of 
0.105(12) e 2 fm 2 derived from this measurement is consistent with those made previously at 60 and 64 
MeV/nucleon, in contrast to an anomalously low result obtained at 43 MeV/nucleon. By coupling a 
multi-configuration description of the projectile structure with realistic reaction theory, the XCDCC 
model provides for the first time a fully quantum mechanical description of Coulomb excitation. 
The XCDCC calculations reveal that the excitation process involves significant contributions from 
nuclear, continuum, and higher-order effects. An analysis of the present and two earlier intermediate 
energy measurements yields a combined B(E1) strength of 0.105(7) e 2 fm 2 . This value is in good 
agreement with the value deduced independently from the lifetime of the l/2~ state in 11 Be, and 
has a comparable precision. 
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The transition strengths between the lowest energy 
states of nuclei are amongst their most fundamental prop- 
erties, being related to the configuration of the states, 
and consequently to particular aspects of the NN inter- 
action. As such, transition strengths are especially im- 
portant for nuclei near the limits of stability, and pro- 
vide important input for nuclear astrophysics. Lifetime 
measurements, whenever possible, are often the preferred 
method for measuring transition strengths. A less direct 
method of determining these strengths is Coulomb ex- 
citation, whereby the projectile is excited via the field 
of virtual photons produced by a heavy target [li]. One 
of the advantages of the Coulomb excitation method is 
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that it can be applied to a wider variety of nuclei. It 
is thus important to validate the method for the well 
known dripline cases. Such work has been performed for 
quadrupole transitions on mostly stable nuclei Q. In 
this paper we focus on dipole transitions applied to an 
unstable system, and resolve a long standing experimen- 
tal anomaly. 

The 11 Be system is an archetype of a halo nucleus 
and exhibits the fastest known dipole transition between 
bound states in nuclei. The B(E1) transition strength 
between the ground and the only bound excited state 
was determined from lifetime measurements by Millener 
et al. to be 0.116(12) e 2 fm 2 Q. 

The anomalously low value for the Coulomb excita- 
tion cross section from a lead target at 43 MeV/nucleon 
4] produced a cross section approximately 40% of that 
expected from Alder and Winther theory of Coulomb 
excitation [5j. Subsequently, further experiments have 
been performed to verify this measurement, utilizing dif- 
ferent targets and energies; MSU at ^60 MeV/nucleon 
on Pb,Au,Be,C [|, RIKEN at 64 MeV/nucleon on Pb 
[7|. The results of these measurements showed broad 
agreement with the lifetime measurement, although at 
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the lower limit of the uncertainties. However, all of the 
subsequent measurements have been performed at signifi- 
cantly (~50%) higher energy than the anomalous GANIL 
measurement. 

Previous analyses of Coulomb excitation experiments 
to extract the 5 (El) of 11 Be USB] have relied upon 
the Alder and Winther theory of Coulomb excitation Q . 
A summary of B(E1) strengths extracted using Alder- 
Winther theory was given in Fig. 2 of Ref. [6j. Although 
this theory gives a useful proportional relationship be- 
tween the -B(El) and the cross section, many approxi- 
mations are applied. The Alder and Winther approach 
assumes that the reaction occurs as a single-step excita- 
tion induced by the virtual photons forming the Coulomb 
field of the target. 

There have been many theoretical efforts to investi- 
gate whether this anomaly could be explained by con- 
tinuum, coupled channel [8| and multi-step effects [9|1, a 
more accurate treatment of the nuclear absorption [lOj. 
and relativistic effects [ll|. As summarized in Ref. [6jj, 
estimates for corrections to Alder- Winther theory at the 
GANIL energy were; continuum effects 4% Q , higher or- 
der effects 6-11% [| and 8% 0], and nuclear effects 8% 
4]. None of these theoretical studies could explain the 
anomalously low cross section of Anne et al. [4| . The the- 
oretical uncertainty from extracting the B(E1) strengths 
using Alder- Winther theory of Coulomb excitation aris- 
ing from all of these corrections is not included in the 
uncertainty quoted (10-13%) 0, d,0|i as this relates to 
the experimental uncertainty only. Additionally, a recent 
analysis of nuclear effects by Hussein et al. [12|, shows 
that estimating the nuclear contribution to the Coulomb 
excitation cross section from the corresponding cross sec- 
tion on the light targets is unreliable. 

The effects of these approximations have been exam- 
ined individually, but a consistent analysis which includes 
all of them simultaneously has so far not been applied. 
One main difficulty is that fully quantum mechanical re- 
action theories, such as the continuum discretized cou- 
pled channels (CDCC) method [HI], treat the projectile 
of interest as a single particle state: in the case of 11 Be, 
a 2si/2 neutron coupled to an inert 10 Be core. Under the 
constraint of a reasonable r.m.s. radius for the ground 
state of 11 Be, a single particle structure model for n Be 
over-estimates the i?(El) transition strength between the 
bound states by a factor of two and thus the predicted 
cross section is too large (see Set Table H]) . Therefore 
in order to extract a reasonable £?(E1) strength from the 
cross section, one must rely on strongly model-dependent 
assumptions. 

Recent advances in reaction modelling [lj| allow for 
an improved structure input in the reaction theory. This 
is the first application of the extended continuum dis- 
cretized coupled channels (XCDCC) method developed 
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in Ref. [14| to Coulomb excitation. The particle-rotor 
structure model, as applied to 11 Be in Ref. yields 
an accurate i?(El) strength which can then be consis- 
tently included in the reaction theory using the formal- 



FIG. 1: (Color online). The Doppler corrected, background 
subtracted 7-ray energy spectrum in coincidence with 11 Be 
ions identified in the 0° telescope, for a n Be beam on a Pb 
target at a mid-target energy of 38.6 MeV/nucleon. The 320 
keV l/2~ to l/2 + transition in 11 Be is clearly identifiable. 
The full line is the result of a Gaussian fit with an exponential 
plus constant background. 



ism developed in Ref. [l4j. Thus a reliable calculation 
for the inelastic cross sections can be made for the first 
time using a fully quantum mechanical description of the 
scattering process. 

Here we report a new experiment performed at 38.6 
MeV/nucleon on a lead target at GANIL. The mea- 
surement presented here was undertaken in an attempt 
to resolve any experimental problems and to elucidate 
any possible energy dependence in the extraction of the 
-B(El). We compare this measurement to theory pre- 
dictions using XCDCC, from which we extract a -B(El) 
strength. We also investigate the effects of various cou- 
plings within the XCDCC framework to explore devia- 
tions from the standard Coulomb excitation theory esti- 
mated in previous theoretical studies [1, [tj [1(1 H3 • Fi- 
nally, we re-analyze results from all previous experiments 
and compare the i?(El) strengths extracted with that ob- 
tained from the lifetime measurement. 

The data presented here were obtained utilizing a sec- 
ondary beam of 11 Be at GANIL, produced by the frag- 
mentation of a primary beam of 18 at 63 MeV/nucleon 
on a 9 Be production target. The fragmentation products 
were purified using the LISE3 spectrometer, resulting in 
a beam of ~5000 11 Be ions per second with 97% purity, at 
a mean energy of 41 MeV/nucleon with an energy spread 
of -1%. The time-of-flight of the ions through the LISE 
spectrometer was measured event-by-event using a PPAC 
located at the entrance to the reaction chamber together 
with an MCP, ~23 m upstream, and also relative to the 
cyclotron RF. These time-of-flight measurements allowed 
the unique selection of 11 Be ions event-by-event from the 
3% impurities in the beam (mainly 9 Li and 13 B). 

The 11 Be beam was incident on a 552 mg/cm 2 Pb 
target, resulting in a mid-target beam energy of 38.6 
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TABLE I: 10 Be-n Woods-S axon potential parameters with en- 
ergies in MeV and lengths in fm, along with the calculated 
£>(E1) strength in e 2 fm 2 , and inelastic cross section in mb. 
Set is the Bel2-pure and Set I is the Bel2-b, both from 
Ref. [l5|]. The experimental cross section is 416(44) mb. 





ffCouiex (mb) 


beam energy (MeV/nucleon) 


39 


59 


Pure Coulomb one-step 


706 


500 


+ nuclear absorption 


579 


421 


+ nuclear in the excitation coupling 


687 


499 


+ bound states coupled to all orders 


642 


469 


+ continuum coupled to all orders 


506 


407 


+ dynamical core excitation (full XCDCC) 


549 


415 



TABLE II: Calculated 11 Be+Pb Coulomb excitation cross 
sections at 39 MeV/nucleon (this work) and at 59 
MeV/nucleon (MSU) showing the effects of the progressive 
improvements in the reaction calculation (see text). 



MeV/nucleon. Charged particles were detected in a 
three-stage telescope centered at 0°, consisting of two 
500 /im thick resistive-strip silicon detectors (providing 
two AE measurements, and measurement of position in 
vertical {x) and horizontal (y) directions), backed by a 
closely packed array of 16 Csl detectors in a 4x4 ar- 
rangement fl6]. The telescope subtended the angular 
range ±9° in x and y. Gamma rays were detected, in 
coincidence with charged particles in the telescope, us- 
ing four Nal detectors at ±45° and ±110° to the beam 
axis, giving a total absolute photopeak efficiency of 3.1% 
for detection of Doppler-shifted 320 keV 7-rays (assum- 
ing an isotropic distribution in the center-of-mass system 
and including relativistic forward-focusing effects of or- 
der 18% and absorption in the target of order 2%). Back- 
ground data were taken with no target present in order 
to determine the contribution from reactions occurring in 
the telescope itself. These were normalized to, and sub- 
tracted from, the target-in data. The present data was 
acquired contiguously with another experiment utilizing 
the same setup [l7j |: the same calibration and analysis 
procedures were used for the present measurement. 

The Doppler corrected 7-ray energy spectrum in co- 
incidence with n Be, after background subtraction, is 
shown in Fig. [T] — a peak corresponding to the 1/2 to 
l/2 + transition (_E 7 = 320keV) is clearly evident, yield- 
ing a cross section of 416(44) mb. The quoted error re- 
flects uncertainties in the fitting procedure, the detector 
efficiencies, the integrated number of n Be beam particles 
and the target thickness. Next we compare with reaction 
theory predictions in order to extract the B(El). 

An accurate prediction of the B(E1) strength can 
be obtained in the particle-rotor model, where n Be is 
treated as a deformed 10 Be coupled to a neutron (see 
Set I Table [T| 15]. This two-body structure model as- 
sumes 10 Be is a rotor and includes its 2 + excited state 
explicitly. It predicts around 15% excited core contri- 
butions in the ground state of 11 Be, in good agreement 
with no-core shell model calculations [18j, and experi- 
ment QJ, ES EH [H. It also produces £(E1)=0.139 
e 2 fm 2 between the two bound states, which is close to 
the value extracted from the lifetime measurement. 

The 11 Be description we use in the XCDCC calcula- 



tions is the model described above. We consider a set 
of parameters based on Bel2-b from Table I of Ref. [i~5| . 
which are given in Set I of Table |TJ Later in the paper 
we will consider the sensitivity of the cross sections to 
various parametcrizations of this potential which yield 
varying _B(E1) strengths. The neutron-target interac- 
tion used is the global optical potential of Schandt et 
al. 12311. The 10 Be-Pb central interaction is taken from 
Ref. [24| . In XCDCC calculations, this interaction is fur- 
ther deformed using /3J — 0.67 for the nuclear part and 
/3§ = 1-13 for the Coulomb, consistent with Ref. [25| . 
When re-analyzing all previous Coulomb excitation data, 
we use the same 10 Be-target potential parameters but 
scale the radius parameters from the Pb to the Au tar- 
get. Errors arising from uncertainties in the nuclear op- 
tical potentials are small [26j . 

Next we give a brief description of the model space for 
the XCDCC calculation. For the n Be bound and contin- 
uum states, we take all partial waves up to £ max = 3. The 
10 Be-n continuum is discretized up to ^p max = 5/2 + and 
E IC \ < 5 MeV. As we are only interested in the inelastic 
cross section to the bound excited state, the details of 
the continuum are not important and thus it is sufficient 
to use rather broad energy bins. For this reason 40 bins 
are included, 5 bins for Jp — 3/2~,5/2 + , with ground 
state core components in the boundary conditions, and 2 
bins for all other partial waves. Radial integrals for the 
continuum bins were performed out to 50 fm with a step 
length of 0.01 fm. 

In Table [TTI we present a summary of the various types 
of reaction calculation that can be performed for this sys- 
tem taking the 11 Be structure model directly from [l5| . 
By turning off couplings within XCDCC, we can provide 
an indication of the size of the various effects which have 
previously been estimated. 

We start with the quantum formulation closest to the 
semi-classical Alder- Winther theory: pure Coulomb one- 
step, taking only Coulomb interactions for incoming and 
outgoing waves, as well as for the transition operator. 
This would be similar to the DWBA formulation with 
Coulomb waves, except that here we go beyond Alder- 
Winther theory as the projectile is described within a 
coupled channels model. The cross section of 706 mb 
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( Table ITU obtained with this simplistic approach is unre- 
alistically large. This is well understood as arising from 
the lack of nuclear absorption: in the Alder- Winther the- 
ory a lower limit on the impact parameter is introduced 
to account for this. The sensitivity of the cross section on 
this is investigated in Ref. [l0(, where a cutoff obtained 
from the particle-rotor structure model is applied. When 
introducing nuclear absorption into the XCDCC calcu- 
lation by including a diagonal nuclear potential in the 
11 Be+Pb optical potential, the cross section is reduced 
to 579 mb. 

The nuclear potential in the distorted waves has a neg- 
ligible effect on the cross section. By subsequently in- 
cluding the nuclear interaction in the transition opera- 
tor, the cross section is increased by 19%. Note that 
nuclear excitation cross section to the 1/2 - state yields 
27 mb. This reaffirms that nuclear effects cannot be sim- 
ply subtracted, as nuclear-Coulomb interference is very 
large [12J. Coupling the two bound states to all orders 
reduces the cross section by 7%. A reduction of 21% is 
obtained with the inclusion of the continuum, resulting 
in a cross section of 506 mb. Finally XCDCC can include 
the dynamical excitation of the core during the reaction, 
which has not been studied previously. This coupling al- 
lows the 10 Be core to change spin through its interaction 
with the target. The ground and excited states of 11 Be 
have different amounts of excited core contributions and 
therefore this coupling allows the relative amounts of ex- 
cited core contributions to re-adjust, enhancing the cross 
section by 9%. Including this effect, which corresponds 
to the full XCDCC calculation, yields a final cross section 
of 549 mb. 

It is remarkable that the full calculation produces a 
cross section similar (within 5%) to the Coulomb-only 
one-step calculation with nuclear distortion (second line 
in Table [TTJ) . This is due to the strong increase arising 
from the nuclear interaction being (coincidentally in the 
current case) nearly canceled by the inclusion of the con- 
tinuum. 

The importance of nuclear effects has little dependence 
on the beam energy and precise scattering angles. An 
analysis of the MSU experiment on the Pb target re- 
veals a similar progression of corrections (+19%,-6%,- 
13%, +2% for an overall change of -2%). The presence 
of effects of the order of 20% suggests again that Alder- 
Winther theory should not be relied upon for the extrac- 
tion of the B (El). 

The calculations presented here are integrated over the 
full angular range. While the MSU measurement had the 
smallest angular acceptance, this still contained 99% of 
the total inelastic cross section. Owing to the significant 
increase in computation time required to get this level 
of accuracy on the angular distribution, the effect of the 
angular cut was neglected for these calculations. For all 
of the other measurements the angular acceptance was 
larger and had a negligible effect on the inelastic cross 
section. 

Relativistic effects are not included in the present 
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FIG. 2: Comparison of -B(El) values obtained from lifetime 
and Coulomb excitation measurements. The weighted aver- 
age of lifetime measurements [3| (open circle) is plotted on 
the left along with the weighted average (solid cicle) of three 
Coulomb excitation measurements (solid symbols). The indi- 
vidial Coulomb excitation measurements, GANIL (this work, 
square), MSU (up triangle) @, RIKEN (down triangle) 0, 
and a previous GANIL experiment (diamond) Q, are plotted 
versus the beam energy. 



work. Relativistic calculations including retardation 
have been performed for 11 Be at 50 and 100 MeV/nucleon 
(llj . Within a semi-classical reaction model, deviations 
of non-relativistic theory from a full relativistic treatment 
were calculated as only 2% at 50 MeV/nucleon, and in- 
creasing up to 14% at 100 MeV/nucleon. Although a full 
relativistic quantum mechanical formulation is still lack- 
ing, for the measurement reported in this paper we there- 
fore expect that errors arising from the non-relativistic 
treatment of the reaction are less than 2%, and around 
5% for the MSU and RIKEN measurements respectively. 

The full XCDCC calculations for 11 Be Coulomb excita- 
tion at an energy of 38.6 MeV/nucleon on a lead target, 
using set I from Table Q] for n- 10 Be, produced a = 549 
mb, somewhat larger than the measured cross section 
of 416(44) mb. As the structure model [15| somewhat 
over-predicts the f?(El) strength (0.139 e 2 fm 2 compared 
to 0.116(12) e 2 fm 2 from lifetime measurements) we ex- 
plore small variations in the parameterization. By vary- 
ing the diffuseness and spin-orbit potential depth, then 
adjusting the odd and even potential depths to repro- 
duce the ground and excited state binding energies, a 
range of fl(El) strengths were obtained (sets I-V of Ta- 
ble UJ . We note that Set III is the parameterization used 
in Ref. [T^. We found a good linear relationship be- 
tween the £?(E1) strength and the Coulomb excitation 
cross section, as would be the case if Alder- Winther the- 
ory were valid. The deviation from the linear regression 
is very small. Finally, we find that the experimental cross 
section is reproduced when the model for 11 Be predicts 
B(E1)=0.105(12) e 2 fm 2 , in good agreement with the life- 
time measurement. 
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We now re-analyze the cross sections from the previ- 
ous experiments and extract B(E1) consistently, using 
the same reaction model, namely XCDCC. The results 
are summarized in Fig.O Our re-analysis of the previous 
GANIL measurement yields the much lower B(E1) value, 
0.048(9) e 2 fm 2 . The authors of this work had mentioned 
in their previous article "this result remains difficult to 
interpret, and even if there is no indication of a possible 
error in the analysis, we feel that it should be regarded 
with some caution until it has been confirmed in a new 
experiment" Q. It has become apparent, based on sub- 
sequent observations of the electronic counting/dividing 
modules employed, that this measurement may have been 
affected by a problem in the normalization of the beam 
intensity [23]. 

For the MSU measurement, we deduce 
S(E1)=0.101(16) e 2 fm 2 and 0.084(12) e 2 fm 2 for 
the Pb and Au targets respectively. For the RIKEN 
measurement, B(E1)= 0.106(11) e 2 fm 2 . All of these 
are consistent with the results from the lifetime mea- 
surement. It may be remarked that the Au target 
measurement, in contrast with the three Pb measure- 
ments, involved a correction of order 50% for 7-ray 
absorption in the target. The three Pb results, with 
quoted errors of 10-15%, are internally consistent at the 
5% level. Taking into account the possible experimental 
problems with the original GANIL measurement and 



also the rather large target correction for the only Au 
data point, the remaining three measurements — all 
with a Pb target — have been combined to give a 
weighted average of B(El)—0. 105(7) e 2 fm 2 . 

Prior to this work, one could speculate that the anoma- 
lous measurement indicated an energy dependence in the 
Coulomb excitation method for the extraction of the 
£>(E1). Our new measurement removes this possibility. 
In summary the work described here demonstrates the 
validity of the Coulomb excitation method for extract- 
ing _B(E1) of exotic nuclei, using a fully quantum me- 
chanical treatment of the reaction process. Moreover, we 
have shown that the accuracy of -B(El) extracted from 
Coulomb excitation can be comparable with that from 
lifetime measurements. 
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